Anderson localization at large disorder 



Laszlo Ujfalusi and Imre Vargcl 

Elmeleti Fizika Tanszek, Fizikai Intezet, Budapesti Muszaki es Gazdasdgtudomdnyi Egyetem, H-1521 Budapest, Hungary 

(Dated: March 4, 2013) 

The localization of one-electron states in the large (but finite) disorder limit is investigated. The 
inverse participation number shows a non-monotonic behavior as a function of energy owing to 
anomalous behavior of few-site localization. The two-site approximation is solved analytically and 
shown to capture the essential features found in numerical simulations on one-, two- and three- 
dimensional systems. Further improvement has been obtained by solving a three-site model. 
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Introduction - The problem of disordered systems has 
been in the front line of research of condensed matter 
physics for several decades starting from the seminal pa- 
per of Anderson^. Electron localization is an essential 
phenomenon playing an important role especially for low- 
dimensional systems^. It is commonly known that the 
spatial extension of the electronic states is a monotonous 
function of disorder and position in the energy-band, 
however, recent studies 3 presented new results that call 
for further understanding. 

In this paper we focus on the large disorder limit of 
the Anderson problem 1 , therefore our Hamiltonian reads 
as 
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where the on-site potential, E{ are uniformly distributed 
over [— -r^, ^] , hence the parameter W characterizes the 
strength of disorder. We measure the energy in units of 
the hopping, t, which is equivalent of setting t = 1. The 
problem described by this Hamiltonian has been stud- 
ied extensively in the past 1 ? 2 . As it is known, there is 
a critical disorder, W c , beyond which, W > W c , every 
eigenstate is exponentially localized. In one and two di- 
mensions W c = 0, meanwhile in three dimensions for the 
above model W c ~ 16.5 for the states in the band center. 

In the present work we investigate the properties of 
the states close to the band-edge in the strongly localized 
regime, W ^> W c . The spatial extent of the eigenstates 
is commonly characterized by the parameter called the 
inverse participation ratio (IPR)£ii, 
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For a state extending homogeneously over k sites I = 
fc _1 , thus 1/1 tells us the effective number of sites a state 
extends to, hence the name. A state localized on one 
single site would give / = 1 , but extending over the whole 
system of size N, I = 1/N. Hence any states will have 
an IPR value between these two cases, 1/N < I < 1. 
In summary the IPR, / is a measure of localization, its 
inverse a measure of extension. 

For strong disorder the states are expected to extend 
over a few sites only, therefore we expect IPR values typ- 



ically of the order of I — >• 1, hence finite size effects will 
not disturb our numerical simulations and in addition rel- 
atively small systems can be used. We employed periodic 
boundary conditions, and linear system size L = 512 in 
ri = 1, L = 20 in d = 2 and L = 8 in d = 3. We computed 
the eigenvalues and the eigenvectors of the Hamiltonian 
(PQ), and the IPR for every state. We made statistical 
averaging over M = 12500 realizations, the results of 
the one dimensional (ID) case are shown in Fig. [1] for 
disorder strength W = 32. In the subsequent part of 
the paper we present our arguments using the ID case 
but extensions for d = 2 and d = 3 will be presented, 
as well. In Fig. QJa) it is shown that the dots are dis- 
tributed symmetrically around the band center, E = 0, 
due to particle- hole symmetry. In Fig. [T] we may also see 
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FIG. 1. (a) The IPR of the states as a function of energy for 
W — 32. The red dots correspond to the states, the black 
curve is the average, (b) The same as (a) but zoomed to the 
band edge. 



that the distribution of the IPRs is broad but many of 
the states seem to have values close to unity, 1=1, and 
there is a considerable amount of them around I = 1/2, 
as well. As a function of energy we can see that moving 
away from the band-center the average IPR, (/) first in- 
creases in accordance with our expectation of increasing 
localization towards the band edge. However, beyond a 
certain energy to be discussed later, 
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the average IPR decreases again, even though the 



2 



Lyapunov-exponent, the inverse of the characteristic 
length of the exponential tail of the state, (not shown 
here) would further increase undisturbed over this en- 
ergy scale. This phenomenon has been recently noticed 
and presented in Ref. y. This means, that at the band 
edge the effective size of a state becomes larger, even if 
the localization length decreases. The main aim of the 
present work is to understand and give an analytical ex- 
planation for the behavior of Fig. [T] especially its part (b), 
showing an empty region of the IPR for large value of the 
energy together with a decrease of the average IPR, (I). 

A qualitative explanation of the existence of this region 
is as follows. If the energy of an eigenstate is bigger than 
W/2 then it must extend at least over two sites, because 
the on-site energies are bounded by \e\ < W/2, hence 
for large energies beyond the potential energy the states 
should have some additional kinetic energy which can be 
obtained by allowing their extension over more than one 
site. This is the reason why close to the band edge the 
states become more extended and hence the IPR has an 
upper bound. A similar argument can be found in 3 . 

In order to understand the large-W behavior of the 
sates first we have to start with the asymptotic behav- 
ior, i.e. as W — >• oo. In that limit all the states will 
be localized to a single site and hopping, i.e. kinetic en- 
ergy, plays no role. In this case the system is a kind 
of sum of independent sites, therefore from the point of 
view of a probabilistic description, it is enough to take 
into account just one site with a random on-site energy 
and one electron on it. The probability distribution for 
E and IPR of this one-site system is the same as for a 
large system due to the independence of sites. The en- 
ergy is just the random potential energy, thus the model 
gives us a W wide band, which is very close to reality 
(see Fig. [1]), for large enough disorder. Since every state 
is localized to one site, for every state 1 = 1 and that 
also becomes a increasingly better approximation as W 
increases. But as mentioned above, for finite disorder 
there is an interesting inner structure in the figures, that 
this one-site model cannot capture. Thus we tried to 
find a better model. For large enough W the states are 
strongly localized and extend over a few sites, which can 
be depicted in Fig. [TJ For example at W = 32 most of 
the IPRs are larger than 1/2, in other words most of the 
states extend approximately over two sites, so a two- or 
three-site model should be enough - at least qualitatively 
- to describe this strongly localized system. 

Models and results - The two-site model As mentioned 
above, an improvement to the asymptotic, large disorder 
limit where the one-site model works, is the so-called two- 
site model. We will show here, that it gives us the main 
physics of Anderson model at large disorder. For such a 
model the Hamiltonian reads as 
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where e± and €2 are uncorrelated random numbers drawn 
with uniform distribution, p(si,£2) = W~ 2 from the in- 



terval [— -y-]. Consequently the support of the prob- 
ability distribution is a square, that is shown in Fig. J^a). 
The eigenvalues and the unnormalized eigenvectors of the 
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FIG. 2. The domain of (a): p(e u £2), (b): p(E ±1 I) for W = 8. 



Hamiltonian ((4]) are 

Si +£ 2 
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Changing s\ and 82 to new variables, t = \{s\ + £2) 
and u = \{e\ — £2), the probability distribution func- 
tion is still constant, but on a square rotated with 45°, 
p(t,u) = 2/W 2 . Using these transformed variables the 
eigenenergies and the IPRs can be writen the following 
way 



£i,2 = t ± \A + u 2 , 
l + 2ii 2 
~ 2 + 2u 2 ' 
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At this point we take the larger eigenvalue, Ei, and ex- 
press t with Ei and u, and then express u using I. After 
these two transformations the probability density func- 
tion of Ei and I become 



p(Ei,I) 
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whose support is quite non trivial: 
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if -\ + 1 < £1 < E and 



< I <l 
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if Eq < Ei < 1 + using the energy border defined 
in Eq. (j3j) and appearing in Fig. [TJ This domain is shown 
in Fig.^b). 
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FIG. 3. p(E, I) (a) The analytical result of the two-site model, 
(b) numerical result for system size L = 512, (c) the differ- 
ence between (a) and (b). (d) Red points are the numerically 
obtained E and I values at system size L = 512, black curve 
is the domain of the two eigenstate in the framework of the 
two-site model. 



The probability distribution function p{E\,I) is ob- 
tained for the larger eigenvalue, E\. For the smaller 
eigenvalue, E2 the result is identical, except E\ must 
be replaced by —E2. Thus the probability density func- 
tion describing the whole system is p(E, I) = p(E\, I) + 
p(E 2 ,I) = p{E\,I) + p(—E\,I), which is shown in 
Fig. [3(a). If we compare this analytical function with 
Fig. [3f b), the probability density function obtained nu- 
merically on a system with linear size L = 512, the qual- 
itative similarity is obvious. In Fig. [I] we can clearly see, 
that the function is symmetric, and it is a sum of two 
components. In our case these two components are the 
two eigenvalues. This simple two-site model also explains 
the behavior of the IPR as a function of energy close to 
the band-edge. Eq. ([TT]) is responsible for the decrease 
of the IPR depicted on the right side of Fig. Jb). An- 
other feature of the two-site model is the peak in p(E, I) 
at / w 1/2. However, there are some differences be- 
tween the model and the numerical results. First of all 
in a two site model / > 1/2, because the state can ex- 
tend maximum to two sites, but in a bigger system there 
exist a few states extending over more than two sites. 
Therefore in the low-/ regime the two-site model natu- 
rally underestimates the reality. Nevertheless both dis- 
tributions are normalized, therefore if somewhere there 
is an underestimation, elsewhere there must be an over- 
estimation, which gives us the hump at high / values in 
Fig. [3(c). Looking at Fig. [3(d) it is clear that the two- 
site model captures very well the shape of the domain of 



definition for the two components, we see a little over- 
estimation for high values of the IPR. From p(E, I) we 
calculated the average of IPR, (/) as a function of the 
energy, E, which can be seen in Fig. [5(b). The analytical 
curve shows a qualitative agreement with the numerical 
function: moving away from the band-center we see an 
increase in (/), and beyond Eq it decreases, showing a lit- 
tle shoulder. Quantitatively in the band-center the model 
overestimates (/), but in the decreasing regime the two 
site model becomes a good approximation (see the inset 
of Fig. 03(b). In addition the two-site model has a band 
of [— ^ — 1, + l] but obviously the real band extends 
beyond these limits. Models and result - the three- site 
model The two-site model introduced in the previous sub- 
section seems to give a qualitatively correct explanation 
for the numerically obtained distributions but as pointed 
out there are deficiencies. In the present subsection we 
will outline the generalization of this model to a three- 
site model. We will investigate how the results change. 
The Hamiltonian incorporating three-sites reads as 
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In the three-site model the domain of the probability den- 
sity function is a cube (see Fig. Ufa)), and the function is 
constant, p(£i,£2>£3) = W~ 3 . As in the two-site model, 
the eigenvectors and IPRs should not depend on the av- 
erage energy, therefore it seems helpful to introduce new 
variables: £\ = t + u, £2 = t + v, 83 = t — u — v. This 
transformation changes the domain to a parallelepiped, 
and the probability density function remains constant, 
because the transformation is linear, p(t,u,v) = 3/W 3 . 
It is easy and straightforward to compute the eigenvalues 
and eigenvectors of Eq. ([T2]) , but the expressions are very 
long, so we do not list here the exact expression, instead 
we only present their support. Every eigenvalue has the 
form Ei = t + Xi( u : v )i with i = 1, 2, 3. Similarly to the 
two-site case, the size of the eigenvectors, i.e. the IPRs 
depend on u and v only, I(u, v). Picking one of the eigen- 
values the problem can be transformed to the variables 
Ei, u and v. The probability density function remains 
constant, p(Ei,u,v) = 3/W 3 . The difficult part of the 
problem is, that the domain changes to a very compli- 
cated object, which is shown in Fig. E(b)(c) and (d) for 
Ei, E2 and E3. To compute p(Ei,I) we have to express 
wasa function of / and v, then calculate the new domain 
and then integrate over v. This resulted in a difficult task 
analytically because the expressions of the IPRs are very 
complicated. Instead we performed our calculation based 
on p{Ei,u,v) = 3/W 3 using a Monte Carlo integration 
over the domains depicted in Fig. 0(b), (c) and (d). The 
result is given in Fig. [5(b). 

In view of Fig. [5(b) it is clear, that the three-site model 
gives a quantitatively better approximation of a large sys- 
tem, especially approaching the edge of the band (see the 
inset), but qualitatively the main behavior is captured al- 
ready by the two-site model. 
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FIG. 4. The domain of definition for variables (a) £i,£2,£3- 
(b) Ei,u,v. (c) E2,u,v. (d) E^^u^v. in a three-site model. 



model and solved analytically capturing the main physics 
of the problem. 

In Ref.0 it is argued that the behavior explained in the 
present work is attributed to the crossover towards reso- 
nant states similar to the effect produced by the Lifshitz- 
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Higher dimensions and summary - Generalization of 
our results to higher dimensions, d = 2, 3 and their com- 
parisons to the numerical simulations are presented in 
Fig. [5] There is a striking similarity between the figures 
which is due to the fact that the strongly localized regime 
is effectively zero dimensional, i.e. as W oo the states 
become localized over a few sites only. The major dif- 
ference is, that the line separating the two components 
becomes less sharp with increasing dimensionality. In one 
dimension this line can be seen very clearly, in d = 2 it 
is still visible, but in d = 3 it becomes hardly visible. 

To summarize we have shown that the Anderson- model 
at strong localization shows interesting behavior espe- 
cially approaching the band-edge. As already known 
the states become more and more localized as energy in- 
creases from the band-center towards the band-edge, i.e. 
the inverse localization length of the states increases as a 
function of energy. The IPR, on the other hand, increases 
up to a critical energy, Eq (j3j) . Beyond this limit the effec- 
tive extension of the states can be described by a multi- 
site (2-site or 3-site) model because in case the eigenergy 
becomes larger than this critical energy, E > Eq, some 
kinetic (hopping) energy is needed besides the random 
potential energy yielding in an upper bound of the IPR 
which in turn results in a decrease of the average IPR, (/) 
as a function of energy in this regime. In order to under- 
stand the numerical simulations we introduced a few-site 
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FIG. 5. Left side: IPR as a function of energy (a) in d = 1 
with L 512 for W 32, (c) in d = 2 with L = 20 for 
W 50, (e) in d = 3 with L = 8 for W = 100. Red dots 
correspond to single states, black curve is the average. Right 
side: (I) as the function of energy (b) in d = 1, (d) in d = 2, 
(e) in d = 3. Dots correspond to the numerically obtained 
curve for a big system, black curve corresponds to the one-site 
model, blue to the two-site model, and green to the three-site 
model. Insets are the same, but zoomed to the left band edge. 



tail^. It would be interesting to find the relation between 
our results and the resonant states. 
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